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ABSTRACT 

A flux splitting scheme is proposed for the general 
non-equilibrium flow equations with an aim at remov- 
ing numerical dissipation of Van-Leer-type flux-vector 
splittings on a contact discontinuity. The scheme ob- 
tained is also recognized as an improved Advection Up- 
wind Splitting Method(AUSM) where a slight numer- 
ical overshoot immediately behind the shock is elimi- 
nated. The proposed scheme has favorable properties: 
high-resolution for contact discontinuities; conservation 
of enthalpy for steady flows; numerical efficiency; appli- 
cability to chemically reacting flows. In fact, for a single 
contact discontinuity, even if it is moving, this scheme 
gives the numerical flux of the exact solution of the Rie- 
mann problem. Various numerical experiments includ- 
ing that of a thermo-chemical nonequilibrium flow were 
performed, which indicate no oscillation and robustness 
of the scheme for shock/expansion waves. A cure for 
carbuncle phenomenon is discussed as well. 

1. INTRODUCTION 

Recently various high-resolution schemes have been 
devised for the Euler/Navier-Stokes equations [1-9]. 
Most of these schemes make use of a first-order up- 
wind differencing as a basis, achieving higher-order ac- 
curacy under some restriction such as Total Variation 
Diminishing(TVD). Since these high-resolution schemes 
inherit their characteristics from their basic schemes, it 
is very important to design a basic scheme with desir- 
able properties. Up to now several basic upwind schemes 
have been proposed, and most of them are categorized 
as either Flux Difference Splitting(FDS) or Flux Vector 
Splitting(FVS). The former uses an approximate solu- 
tion of the local Riemann problem, while the latter splits 
the flux vector into up-stream and down-stream travel- 
ing components. 

The most rigorous approach of the FDS scheme may 
be the use of the exact Riemann solver, as proposed 


by Godunov[10]. Since the Riemann problem is a non- 
linear problem, this scheme requires iterative numeri- 
cal procedures, leading to an inefficient method espe- 
cially for vector machines or for gases with a general 
equation of state. Because all that is needed from the 
Riemann solution is a flux across a cell interface, some 
approximations should be possible. The most success- 
ful and popular approximate Riemann solver is Roe’s 
method, which is the exact solution of a linearized Rie- 
mann problem[ll]. This is able to capture any single 
stationary discontinuity with no numerical dissipation, 
unfortunately including an expansion shock. Such an 
unphysical solution can be removed by “ Entropy Fix[5]” 
At strong expansions, however, the Roe scheme diverges 
even if the entropy fix is applied. In addition, it is 
known that the Roe scheme has a more serious prob- 
lem, so-called “carbuncle phenomenon” which is a nu- 
merical instability in capturing a strong shock wave 
in multi-dimensional computation. Recently Quirk [12] 
has reported that even the exact Riemann solver suf- 
fers from this problem. Roe’s approximate Riemann 
solver has consistency with the exact one about the 
conservation law, and this type of approximate Rie- 
mann solvers are called “ Godunov- type” schemes [13]. 
Another Godunov-type scheme is the HLLE scheme, 
which much simplifies the Riemann problem by express- 
ing the solution with only two waves. Einfeldt[14, 15] 
has shown that this scheme satisfies the positivity and 
entropy conditions if some appropriate velocities are 
chosen for these two waves. The HLLE scheme, how- 
ever, lacks the information of contact discontinuity, so 
it is too dissipative to be applied to the Navier-Stokes 
equations although that dissipation is strong enough 
to remove the carbuncle phenomenon[12]. The HLLE 
scheme can be modified to include the information of the 
contact discontinuity under the positivity condition [16]. 
However, that scheme again suffers from the carbuncle 
phenomenon [17]. Also numerical noise radiated from 
slowly-moving shock waves has been reported for the 


1 



exact Riemann solver as well as the Roe scheme[18], 
and is probably a common problem for any Godunov- 
type schemes. Another approximate Riemann solver 
was devised by Osher[19]. This scheme is not rigor- 
ously a Godunov-type scheme, but it captures a single 
shock wave at most by two-intermediates points, and 
reproduces the exact solution if a single expansion wave 
or contact discontinuity is present. Further, Osher’s 
scheme does not produce the noise at slowly-moving 
shock waves. However, it fails in near vacuum condi- 
tions because it uses the reversed-order integration path 
in phase space, which is needed to avoid a serious over- 
shoot in capturing shock waves. 

In FVS, the flux function F is divided into positive 
and negative parts: 

F(q) = F(q)- + F(q) + , (1) 

which give the numerical flux at cell interface between 
the states qi, and qj*: 


convective quantities. The AUSM can capture a station- 
ary contact discontinuity with no numerical dissipation 
and is robust enough to calculate strong shock waves. 
However, it bears a slight numerical overshoot imme- 
diately behind the shock. Roe’s assertion in 1988 [28], 
“The perfect Riemann solver (including FVS) for gener- 
ating numerical schemes does not yet exist” seems still 
true. 

In this paper, we studied another way to remove the 
numerical dissipation of the Van-Leer-type flux-vector 
splittings on a contact discontinuity. The scheme ob- 
tained is also recognized as an improved AUSM scheme. 
Various numerical experiments, including a thermo- 
chemical nonequilibrium flow problem, are performed. 
A cure for the carbuncle phenomenon is discussed as 
well. Throughout this study we basically concentrate 
on the explicit first-order scheme in space and time in 
order to remove any uncertain factors other than the 
numerical flux itself. 


*\/ 2 (qi,i q*) = F(qi f ) + + F(qj*) (2) 

Steger and Warming [20] made use of similarity trans- 
formations and the homogeneity property of the Elder 
equations, splitting the flux depending on the sign of 
eigenvalues of the flux Jacobian matrix. This splitting 
showed errors around the sonic point, and Van Leer [21] 
proposed an alternative splitting, which gives noticeably 
better results and produces steady shock profiles with at 
most two zones. Hanel[22] further modified Van Leer’s 
FVS so as to preserve total enthalpy for the steady state 
solution. These FVS schemes, however, share a serious 
disadvantage, that is, large dissipation on contact dis- 
continuities and shear layers. Liou and Steffen [2 3] tried 
to remove the dissipation by using higher-order polyno- 
mial expansions of the flux vector, but their scheme was 
not robust enough in multi-dimensional calculations. 

In another effort to develop less-dissipative upwind 
schemes, the surplus dissipation of the FVS is reduced 
by introducing the flavor of FDS into FVS schemes. 
These schemes are not FVS any more because their nu- 
merical flux can not be expressed in a form of Eq.(2), 
and we call them FV /DS. Hanel[24] has found that the 
numerical dissipation in the boundary layer is greatly 
reduced by using one-sided upwinding of the tangen- 
tial velocity in Van Leer’s FVS formulation. This w as 
later extended by Van Leer[25], who employed one-sided 
upwinding also for the enthalpy of the energy flux. But 
their schemes still hold the numerical dissipation for 1-D 
contact discontinuities and yield glitches in the pressure 
near the edge of the boundary layer. Liou and Steffen 
proposed a more promising scheme named Advection 
Upstream Splitting Method(AUSM)[26, 27], in which 
the cell-face advection Mach number is appropriately 
defined to determine the upwind extrapolation for the 


2. GOVERNING EQUATIONS 

2.1 Generalized Nonequilibrium Flow Equations 
Generally, a nonequilibrium flow such as described 
by chemical reactions or turbulent models has the gov- 
erning equations in a form: 


dq dF k _ 
dt dx k 


( 3 ) 
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and 

J) = p(/>,C,p/i,p/2,...,P/n)- ( 4 ) 


These equations include the conservation of total mass, 
momentum, total energy, and also the physical quanti- 
ties p/*, which represent nonequilibrium effects, t.e., ei- 
ther the concentration of chemical species or vibrational 
energy. The quantities p, p, u k , E , and e respectively 
denote the pressure, density, Cartesian velocity compo- 
nents, total and internal energies, whereas the vector S 
is a set of elements of nonequilibrium source terms. It is 


2 


assumed that an appropriate “/rozen” speed of sound c 
is calculated depending on a gas model, by which Eq.(4) 
is defined. In this study, a numerical scheme is formu- 
lated for the governing equations given by Eqs.(3) and 
Eq.(4). 

2.2 Two- Temperature model 

In this work, Park’s two-temperature model 
equations[29] are used for calculating nonequilibrium 
flows. This is among one of the simplest physical 
models for describing the thermal-chemical nonequilib- 
rium phenomena. It employs an approximation that 
the translational/rotational energy modes are in equi- 
librium at the translational temperature, while the 
vibrational/electronic/electron-translational modes are 
at the vibrational temperature. The conservative vari- 
ables of this model include the mixture of vibrational- 
electronic energy, ey , as well as those of concentration 
of chemical species: 

q = (p,pui,pU2,pU3,E,pzi,pZ2,...,pz n ,evY, (5) 


where z s is molar number density per unit mass of 
species s. The chemical reaction rates depend not only 
on the translational temperature but also on the vibra- 
tional temperature, Ty , which is related to the mixture 
vibrational-electronic energy ey: 


ey = ^2 pz, e V ib , (Ty ) + x jTpz, e e u , (Ty )+pz e - RTy . 


szzmol 


( 6 ) 

The vibrational energy e v «6« per unit mole is obtained 
by assuming a harmonic oscillator, and the electronic 
energy e e j c * is calculated by taking accounts of only the 
first electronic excited state. The conservation equation 
of ey is given as 


dey dey u* 

sr + -sir = nr ’ 


( 7 ) 


where 


sy 


Y] + 

a —moi 
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g—mcl 


e vib» ~ 
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du k 




h ® [( w+ «,)g-l + 


dr , 


The first term in sy is the vibrational energy lost 
or gained due to molecular dissociation or combina- 
tion. Here, w $ is chemical production rates of molec- 
ular species s, and D t is average vibrational energy per 
unit mass of molecule s, which is created or destroyed 
at rate w t . The second term is the energy exchange be- 
tween vibrational and translational(T-V) modes due to 


collisions with el ihs being the vibrational energy in equi- 
librium with the translational energy. Usually an ex- 
perimental formula in [30] is used as relaxation time r, , 
which, however, has been found to predict a too fast T- 
V relaxation at high temperature. Park introduced two 
modifications to remedy it. One is the collision limiting 
of the vibrational relaxation times, and the another is 
the diffusive T-V relaxation [29]. Another formula which 
takes accounts of high-temperature effects on the T-V 
relaxation is found in [31]. The third term in sy is the 
work on electrons done by an electric field induced by 
the electron pressure gradient. The forth term is the en- 
ergy exchange due to elastic collisions between electrons 
and heavy particles where u t9 is collision frequency for 
electrons and heavy particles. The remainder terms are 
due to the conduction and diffusion. More detailed de- 
scription of the two-temperature model is found in [32]. 
In this study, the 7-species chemical reaction model[33] 
was used: 

(1) 0 2 +M^O + 0 + M, M=N, NO, O, 0 2 , N 2 ; 

(2) Na + M^N + N + M, M=N, NO, O, 0 2 , N 2 ; 

(3) NO + M- N + 0 + M, M=N, NO, O, 0 2 , N 2 ; 

(4) 0 + N0^N + 0 2 ; 

(5) O + N 2 ** N + NO; 

( 6 ) O + N ^ NO+ 4- c“. 

2.3 Numerical Flux in Generalized Coordinates 
In practical computation, a numerical flux in gen- 
eralized coordinates is needed, which we calculate af- 
ter [7]. Let the vector n = (ni,n 2 ,n 3 ) be a normal- 
ized cell-interface normal in ^-direction with the vectors 
and m(mi,m 2 ,m 3 ) being its recipropal ones: 
n • 1 = 0; n • m = 0; 1 • m = 0; |n| = |1| = |m| = 1. The 
normal and tangential velocity components to the cell- 
interface normal are calculated for each left and right 
state: 

u = niUi, v = kui i w = miUi. ( 8 ) 

The numerical flux in the x-direction of this local Carte- 
sian coordinates has the form: 


pu 

pu 2 +p 
pvu 
pwu 
pHu 

pfi* 

pf 2 u 


pfttU J 


where H is the total enthalpy: 


ff = £±f. 

P 


( 9 ) 


( 10 ) 


3 



Following a specified flax splitting scheme, the numeri- 
cal flux Fx/j = (Fi, Fa, . . . , Fs+»)‘ is calculated in these 
local Cartesian coordinates. Finally the numerical flux 
in the ^-direction of the generalised coordinates, F{, is 
given as 


F e = 5 


Fi 

ni Fa + Jifj “I" m if< 

naFa + fa-fa 4- 1 H 3 F 4 

TI 3 F 3 + I 3 F 3 + TH 3 F 4 

F 5 

F e 


( 11 ) 


L fs+n J 

where S is the area of the cell interface. In this proce- 
dure the numerical flux in the generalized coordinates is 
uniquely specified by the definition of the numerical flux 
in the local Cartesian coordinates, which is discussed in 
the next section. 


3. NUMERICAL SCHEME 

3.1 Classification of the PV/DS schemes 

Up to now, the most successful FV/DS schemes have 
been the Van Leer scheme [2 5] and the AUSM[26, 27]. 
These schemes are equipped with favorable properties: 
conservation of enthalpy for steady flows; small dissipa- 
tion in the shear layer. As the first step in constructing 
our scheme, we generalize these schemes, calling them 
A US MV- type and AUSMD-type schemes, where “V” 
and tf D w denote flux- Vector-splitting-biased schemes 
and flux-Difference-splitting-biased ones, respectively. 

It is useful to first review briefly the AUSM scheme. 
An interesting form for writing the numerical flux is 
given in [27], as 

Fi/a = x[ + 1 + Pi/2> 

( 12 ) 

where 


# = {p,pU,pV,pW,pH,pfl,... t pfnY, 

Pi/2 = (0. Pi/2| o, 0, 0, 0, ... , 0)‘. 

Here the interface velocity ux/? and pressure pi /2 are 


defined as 
where 


u* = 


and 


ttl/2 = «£ + 

(13) 

±^-(«±c) J , if>[<c; 

1 

-(u ± |u|), otherwise, 

z 

(14) 

P 1/2 = Pl + Pr< 

(15) 


where 

_± _ / p(“/c ± 1) 2 ( 2 T «/<0/4> if M < c; , 

f \ p(u± |tt|)/2tt, otherwise. ' * 

Also Liou and Steffen previously proposed another 
AUSM scheme, which is based on Mach number 
splitting[26], where an appropriately defined interface 
Mach number is used to determine the upwind extrap- 
olation for the convective quantities. 

Now, we define the AUSMD-type schemes, which in- 
clude both of the AUSMs as a member. In defining this 
type of schemes, we use the interface mass flux, {pu)\/ 2 , 
instead of the interface velocity or Mach number as a 
switching function for the upwind extrapolation. The 
AUSMD-type schemes are defined by their numerical 
flux: 

F 1/2 = (pw)i/2(*i+*«) - l(p«)i/al('®Ji-*i) ]+Pi/a> 

(17) 

where 

¥ - (1, u, r, w y H t / 1 , . . . , / n )‘. 

A specific numerical scheme is uniquely defined by an 
appropriate interface mass flux {pu)\f 2 and interface 
pressure P1/2 • So far, the interface pressure given by 
Eq.(15) and Eq.(16) has been commonly used in the ex- 
isting FV/DS schemes, but the form of the mass flux 
varies with each numerical scheme. This is because 
the mass flux is directly connected to the resolution of 
contact discontinuities, which is a main interest in the 
FV/DS schemes. For example, the following mass fluxes 
may be possible for the AUSMD-type scheme, and are 
actually used in Van Leer’s FV/DS and the AUSMs, 
respectively. 

• Van Leer’s FV/DS: 

(pu)i/ 2 = p L + p R . (18) 


• AUSM (velocity-splitting): 

(pu)i/2 = Uii2{pL +Ar) “ l^i/al {pr-Pl)), 

(19) 

where 

u l/2 = U L + U R* 

• AUSM(Mach number-splitting): 

(pu) 1/2 = M i M plcl + PrCr) 

- |Afi/ 2 | {prcr - Plcl) ]■ (20) 

where 

■W 1/2 = + Vr/ c R- 
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Here are defined by Eq.(14). The interface mass flux 
defined by Eq.(19) or Eq.(20) in conjunction with the 
interface pressure of Eq.(15) makes the AUSMD-type 
scheme reduce to the velocity-splitting-based AUSM 
and the Mach-number-splitting- based AUSM, respec- 
tively. It is noted that the mass flux of the Mach- 
number-splitting-based AUSM vanishes at a station- 
ary discontinuity, making the AUSM a less dissipative 
scheme among the FV/DS schemes. In the next section, 
we will present a new formula for the interface mass flux 
(pu) 1/2 which results in noticeable improvements over 
the AUSM. 

The numerical flux of the Van Leer FV/DS 
scheme [2 5] slightly differs from the form of Eq.(17) even 
if the mass flux of Eq.(18) is used, because that scheme 
does not use the upwind extrapolation about the term 
(pv?)i /2 in the x-momentum flux. Hence, we need to de- 
fine another class of FV/DS schemes — AUSMV-type 
schemes: 

(P^AUSMV = u t + U R (*"*)*• ( 21 ) 

But the AUSMD-type schemes defines 

(P^AUSMD — ^[(/"Oi/afalr+i 1 *) - |(pu)i/ 2 |(uii-ui,)]. 

( 22 ) 

The mass flux of Eq.(18) and the velocity/pressure split- 
tings of Eqs.(14,16) make the AUSMV-type scheme re- 
duce to the Van Leer FV/DS scheme. We emphasize 
that a specific AUSMV-type scheme needs the defini- 
tion of the velocity splitting, u ± , for the x-momentum 
flux as well as that of the interface mass flux and the 
pressure. The characteristics of the AUSMD-type and 
AUSMV-type schemes will be discussed in detail in the 
following sections. 

3.2 Removal of Numerical Dissipation at Contact 
Discontinuities 

We discuss first the formula for the interface mass 
flux because both AUSMD-type and AUSMV-type 
schemes need it. The main drawback of Van Leer’s 
FV/DS scheme is the numerical viscosity on the contact 
surface, while that of the AUSM is the numerical over- 
shoot at shock waves. We have found that the overshoot 
of the AUSM is mainly due to the mass flux differenc- 
ing. This is probably because the AUSM mass flux of 
Eq.(19) nor Eq.(20) does not directly take into account 
of the density behind the shock wave. Hence, in this 
study we employ a mass flux formula of Eq.(18). In this 
case, however, the surplus numerical dissipation at the 
contact discontinuities will be a problem. 

First, let’s consider a stationary contact discontinu- 
ity where ul = ur = 0 ;px, = PR • The mass flux of 
Eq.(18) is 

{pu)if2 = (plcl - Prcr)/*> ( 23 ) 


which does not necessarily vanish, producing the exces- 
sive numerical dissipation at the contact discontinuity. 
The term of ci and cr in Eq.(23) comes from the numer- 
ical dissipation which is added other than the simplest 
velocity splitting: u ^ = (u ± |u[)/2. Some additional 
dissipation, however, is indispensable in order to get a 
stable scheme. Hence, we modify the velocity splitting 
so that the numerical dissipation can cancel at the con- 
tact discontinuity: 


4. r (« + c) 2 t*+M, , «i + l«i|. 

*i =ocl{ — c ~ h+ 2 / 


= «r{— 


(« - c) 2 U - |u|-, , Ur - |«fl| 


where 


and 


<*i = 


4c 


2 & L 


&L+&R 


<XR 


'-}r + 


2&r 

&L+&R 


<*L = 


/ (PL i PR ) 

(pc)i ’ 


OR = 


/(Pfl.Pi) 


(24) 

(25) 

(26) 


and “/” is some function of pressure pl and pr- Since 
p L equals to pr at the contact discontinuity, an arbi- 
trary function might be possible as /. According to our 
numerical experiments, we take 


/(Pi.P*) = Pi. /(P*.Pi)=P * ■ ( 27 ) 


The velocity splitting, which is a function of a: 0 < a < 
2, is shown in Fig. 1. 
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This velocity splitting and the interface pressure de- 
fined by Eqs.(15,16) make both of the AUSMD/V-type 
schemes satisfy the condition of no dissipation at the 
contact discontinuity on every component of the numer- 
ical flux vector. It is noted that, if the velocity splitting 
of Eq.(24) is applied to Van Leer’s FVS[21], that scheme 
is also able to keep the stationary contact discontinuity. 

Now, let’s consider a moving contact discontinuity: 
ul = ur = t ilr 0, where the exact solution of the 

Riemann problem gives a mass flux: 


(pu)i/3 = vlr(pl + PR.) - |«X,Jl|(pfl - Pl) ]• (28) 

In this case, the use of a common speed of sound, c™, 
instead of cj, and cr in the modified velocity splitting 
of Eq.(24) enables a mass flux formula of Eq.(18) to 
reproduce that of Eq.(28). The condition on the mo- 
mentum flux is also satisfied by the use of the com- 
mon speed of sound in the pressure splitting of Eq.(16). 
This simple modification makes all other fluxes repro- 
duce those of the exact Riemann solution for the con- 
tact discontinuities. In this study, max(c£,cji) is taken 
as the common speed of sound. We call the spe- 
cific AUSMD/V-type scheme which has the above in- 
terface mass flux and pressure/velocity-splittings as the 
AUSMD and AUSMV, respectively. 

We now summarise the numerical scheme for the 
AUSMD/V. First, the mass flux is 

(ptt)i/ 2 = « t PL + PR> ( 29 ) 

where the velocity splittings , ti* are no longer the 
familiar Van Leer splittings, but rather including terms 
designed to capture stationary/moving contact discon- 
tinuities, and are given as 


f a L { ifl2 £l<l; 

I m LtbiLL, otherwise, 

(30) 

uZ = 


i 


tt *~h**l otherwise, 

2 1 

(31) 


where 

= 

and 


2 ( p/p)l 


( v/p)l + (p/ p)r ’ R (p/ p)l + (p/ p)r ' 


2 (p/p)r 


Cm = max(ci.Cft). 

Secondly, the pressure flux is 

Pi/j = Pl +Pr> 


(32) 

(33) 

(34) 


where 


[ riMm? ± i)‘p t ^ < i; 

»!/«= otherwise. 

I 2 ttj ,/R 

(35) 

Substitution of Eq.(29) and Eq.(34) into Eq.(17) results 
in the numerical flux of the AUSMD, in addition, if the 
term pv? in the x-momentum flux is replaced by Eq.(21), 
we have the scheme AUSMV. 

3.3 AUSMD V: Mixture of AUSMD and AUSMV 

There remains a question about the choice between 
the AUSMD and AUSMV, t.e., between Eq.(22) and 
Eq.(21). This will be decided by numerical experiments. 
First the shock capturing capability is tested by solving 
a shock-tube problem for ideal gas with specific heat 
ratio, 7, of 1.4. The initial condition is (p,p, M)l = 
(0.1,0.1,15), (p, p, M)r - (0.1, 0.1, -15), which pro- 
duces strong shock waves propagating in both direc- 
tions. Figures 2(a), 2(b) and 2(c) show the profiles of 
temperature, pressure and density after 100 time steps 
at CFL=1.0 with 100 grid points by the AUSM, AUSMD 
and AUSMV, respectively. The AUSM and AUSMD 
show noticeable oscillations, although the AUSMD is a 
little better than the AUSM. On the other hand, the 
AUSMV gives a fairly good solution in Fig. 2(c), indi- 
cating that the AUSMV has a higher shock-capturing 
capability than the AUSMD. 

Another shock-tube problem is a slowly moving 
contact discontinuity: (p»Pi*i)/e = (10.0 , 1.0, 0.3cr), 

(p, p, u)l = (0.125, 1.0,0.3c*) with 30 grid points. Fig- 
ure 3(a) shows the solution after 100 time steps by the 
AUSMD. Figure 3(b) is the solution by the AUSMV. 
The AUSMV produces serious oscillations, which take 
place only on the velocity (momentum) field a inside ” 
the left state but not at the contact discontinuity. To 
clarify the cause of the oscillations, a simplified prob- 
lem is set up, in which a center grid point is given 
a small velocity perturbation (u = 0.001c) from a uni- 
form field: (p,p, u) L / R = (1, 1, 0). Figures 4(a) and 4(b) 
are the results after 15 time steps with CFL number of 
0.75 and grid number of 31 by the use of the AUSMD 
and the AUSMV, respectively. For this simple prob- 
lem, AUSMV produces serious oscillations, although it 
works well when CFL number is reduced to 0.5. This 
kind of oscillations, however, are a common problem 
for the schemes which use a momentum flux formula of 
Eq.(21): Van Leer’s FVS; Hand’s FVS; Hand’s FV/DS; 
Van Leer’s FV/DS. The vdocity at the perturbed grid 
point at the second time step is analytically derived for 
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various schemes, and is 

(1 — A)c :Roe’s scheme; 

(1 _ i_A)e :AUSMD, AUSM; 

2y 

(1_HL±1a)c :AUSMV, 

Van Leer/Hanels’ FV(/D)S, 

(36) 

where A = cAt/Ax is the CFL number. This analysis 
indicates that a small CFL number< 27 /( 27 - 1 - 3 ) ~ 0.48 
should be used for the AUSMV scheme in order to avoid 
the overshoot in the velocity field. This restriction seems 
too strong, even though the AUSMV scheme shows a 
nice shock capturing capability. Hence, we prefer a 
mixed momentum flux of (pu)ij 2 between the AUSMV 
and AUSMD: 




Fig. 2(b) Shock-tube problem — colliding flow by AUSMD. 


(pv 2 ) 1/2 = d + *)(p« 2 )ausmv + (5 - »)(^ j )ausmd. 

(37) 

where $ is a switching function of the pressure gradient: 


1 . ^ I Pr-Pl\ ^ 

s — - min(l, iv r). 

2 mm(pi,,PjO 


This averaging is biased toward to the AUSMV in order 
to maintain the shock-capturing capability. We call this 
mixed scheme the AUSMD V. In this study a constant 
parameter K = 10 is taken. The AUSMD V is applied to 
the previous shock-tube problem for colliding flow. The 
result is shown in Fig. 2(d), which indicates that the 
AUSMDV has a robust shock-capturing capability. It is 
noted here that the genuine AUSMV flux is also possi- 
ble if the CFL restriction is not considered serious. The 
results for the moving contact discontinuity problem by 
the use of the AUSMDV and other schemes are found 
in the Section 4.1. 


AUSWM-Splitting) CFL-1.00 N-100 



Fig. 2(a) Shock-tube problem — colliding flow by AUSM. 


AUSMV :APuP^Pt CFL-1.00 N-100 



Fig. 2(c) Shock-tube problem — colliding flow by AUSMV. 



Fig. 2(d) Shock-tube problem — colliding flow by AUSMDV . 


7 





Fig. 3(a) Shock- tube problem — moving contact disconti- 
nuity by AUSMD. 



Fig. 3(b) Shock-tube problem — moving contact disconti- 
nuity by AUSMV, 


AUSMD :APuP*)-Pl 


CFU0.75 N- 15 

— * — Pressure 
— w — Density 
• 'Temperature 
— — u x6 + 1 


0 3 


l 1 > 0.0 - 

. 1.0 - 0.5 0 0 



Fig. 4(a) Shock-tube problem — small perturbation on ve- 
locity field by AUSMD. 


Fig. 4(b) Shock-tube problem — small perturbation on ve- 
locity field by AUSMV. 


3.4 Shock Fix — a Cure for the Carbuncle Phe- 
nomenon 

Recently, Quirk[12] has proposed a test problem to 
investigate the carbuncle phenomenon. The test prob- 
lem is a shock wave propagating into a static gas through 
a duct whose centerline grid is slightly perturbed: 


f Y jtk mid = Yhmid + £ for j even, 

\ Y jtkm id = Y k mid - e for j odd. 


In this study, a unit-spaced 350x19 grid system is used 
with the central grid perturbation t being 10~* 6 , and the 
shock Mach number is 6. 

Figure 5(a) shows the density contours solved by 
Roe’s approximate Riemann solver after 400 time steps 
at CFL number of 1.0, in which the typical carbuncle 
phenomenon takes place. The same problem was solved 
by the AUSMDV. That result is shown in Figure 5(b), 
where the carbuncle phenomenon also emerges although 
it is less serious than that of the Roe scheme. Other 
schemes: the exact Riemann solver; the Osher scheme; 
the HLLE scheme; Van Leer’s FVS; Hand’s FVS; the 
Hand/Van Leer FV/DS; the AUSM, were tested for 
this problem. Surprisingly, most of these schemes: the 
exact Riemann solver; Osher ’s scheme; the Hanel/Van 
Leer FV/DS, suffer from the carbunde phenomenon, al- 
though degree of instability may vary. It is noticeable 
that the AUSM is free from the carbuncle phenomenon 
even though it is a less dissipative scheme. However, 
the carbuncle phenomenon seems to be a rather com- 
mon problem for usual shock- capturing schemes. 
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Fig. 5(a) Shock wave propagating through a duct by Roe 
scheme. 



Fig. 5(b) Shock wave propagating through a duct by AUS- 
MDV. 

We suspect that the carbuncle phenomenon comes 
from the recognition of a multi-dimensional shock wave 
by a one-dimensional numerical scheme. When the 
shock wave is captured by a shock capturing scheme, 
an intermediate point (s), which is unphysical, is needed 
to express a shock position numerically. In 1-D case, if 
the scheme is well designed, this is not a problem. But 
when that scheme is applied to multi-dimensional cal- 
culation, the intermediate point exchanges information 
with the neighbors which are also intermediate shock 
points. It is probable that the exchange of information 
between these unphysical data causes a numerical insta- 
bility. In that case, a mechanism to stabilize the shock 
or some form of artificial viscosity is required. We call 


this artificial procedure “Shock- Fix.” Here, the follow- 
ing shock-fix is proposed. 

1. Find out compressible sonic points across the cell 
interface before calculating numerical flux. This 
sonic point takes place when 

{ ul-cl > 0 and u R - c R < 0, 
or uj, + cl > 0 and u R + c R < 0. 

If this is detected between grid point j and j + 1 , 
set flag S^j = 1 and S^j+i = 1 in ^-direction. The 
same procedure is done in T)/(- direction, which is 
easy to vectorize. 

2. In calculating a numerical flux in ^-direction be- 
tween grid points j and j + 1, S v and but 
are used to sense the shock position, to which some 
dissipative scheme is applied. That is 

Dissipative Scheme: if 

(S v ,j + S Vi j+i + S( t j + 1 ) > 1 ; 

Non- Dissipative Scheme: otherwise. 

The numerical fluxes in other directions are calcu- 
lated in a similar way. 

This method unlike [12] does not use the pressure gradi- 
ent in sensing shock waves. The pressure gradient does 
not always provide a sufficient information for the recog- 
nition of shock wave, since in that case some empirical 
parameter to distinguish shock and compressive wave 
would be needed. It is noted that the alone shock fix 
procedure is easily applicable to unstructured grid ap- 
proach. This shock-fix needs a dissipative scheme to 
stabilize the shock wave. For the AUSMDV, we employ 
Hand’s FVS scheme, because it is dissipative enough 
and conserves the total enthalpy for steady flows. The 
numerical flux of Hand’s FVS is 

Fi/2 = PL^t + PRUji^R + Pl/2> ( 40 ) 

where the velocity splitting and the interface pressure 
are defined in Eq.(14) and Eqs.(15,16), respectively. 

Figure 5(c) shows the solution by the AUSMDV 
with the proposed shock-fix, where the carbuncle phe- 
nomenon is removed now. A similar shock-fix is pos- 
sible for Roe’s approximate Riemann solver, for which 
the HLLE scheme is considered to be a dissipative part- 
ner scheme. The solution of the Roe scheme with this 
shock-fix is shown in Figure 5(d), where the instability 
is totally eliminated. This indicates the validity of the 
proposed cure of the carbunde phenomenon, and seems 
to support our suspicion about its origin. 
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Y-axis( behind shock) X-axis(centerUne) 


Fig. 5(c) Shock wave propagating through a duct by AUS- 
MDV with shock-fix. 



Fig. 5(d) Shock wave propagating through a duct by Roe 
scheme with shock-fix. 

8*5 Entropy Fix 

It is well known that the Roe scheme without 
entropy-fix gives an expansion shock. But also the Go- 
dunov scheme(exact Riemann solver) produces glitches 
at sonic points [34, 35]. To investigate this problem, a 
shock-tube problem with an initial condition of the in- 
verse shock of Mach number 15 was solved. Figures 
6(a) and 6(b) show the density profile on a mesh of 
100 points after 50 time steps at CFL number of 1.0, 
using the Roe scheme and the Godunov scheme, respec- 
tively. The Roe scheme keeps this unphysical disconti- 
nuity, while the Godunov scheme spreads it out but still 
has a slight glitch at the sonic point. The AUSMDV 
scheme also produces a similar glitch in Fig. 6(c). 


Ro« CPU 1.00 N- 50 



Fig. 6(a) Shock-tube problem — inverse shock wave by Roe 
scheme. 


Exact Riemann Sohrer C FI- 1.00 N- 50 



Fig. 6(b) Shock-tube problem — inverse shock wave by Go- 
dunov scheme. 

AUSMDV : Apuf^p^lO CFL- 1.00 N- 50 



Fig. 6(c) Shock- tube problem — inverse shock wave by AUS- 
MDV. 

Usually the numerical flux near the sonic point needs 
some modification to remove this glitch. In this study, 
the following modification is proposed. There are two 
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cases where an expansive sonic point is found across the 
cell interface: 

f Case A: u L - c L < 0 and ur - cr > 0; 

{ Case B: ujj + cx, < 0 and uj* + cr > 0. 

Then we introduce a numerical dissipation only when 
a single expansion wave is detected, modifying the flux 

F 1/2 to Fi/ 2 
if( A and B ) then 

no modification 

elseif ( A ) then 

P = Pl/2 - CA(« - c)A(p¥) 
elseif ( B ) then 

= ®*i /2 — CA(ti + c)A(p®) 

else 

no modification 

endif 

where A( ) = ()*-( )l, and a constant parameter 
C is taken 0.125. This entropy-fix keeps the enthalpy- 
conservation property of the AUSMDV. Figure 6(d) 
is the shock-tube solution by the AUSMDV with this 
entropy-fix, where the glitch is fairly removed. 


£U$MDV :ftj>uPJ=PLJ=*0 CFL-1.00 N- 50 



Fig. 6(d) Shock-tube problem — inverse shock wave by AUS- 
MDV with entropy fix. 

Another numerical example is a strong shock 
diffracting around a corner[12]. Figure 7(a) and 7(b) 
show the density contours by the Godunov scheme and 
the AUSMDV without entropy-fix, respectively. These 
solutions are similar, and give rise to a nearly discontin- 
uous expansion fan. Figure 7(c) is the solution by the 
AUSMDV with the entropy-fix, where the discontinuity 
is removed fairly well. The inner shock wave of the AUS- 
MDV looks weaker than that of the Godunov scheme, 
but its strength seems to be related to the glitch of the 
expansion fan, for the Roe scheme with the entropy- 
fix[34] gives almost the same solution as the entropy- 
fixed AUSMDV in Fig.7(c). It is noted that the Roe 
scheme without the entropy-fix and the Osher scheme 
failed in this problem. 



Fig. 7(a) A shock wave diffracting around a comer by Go- 
dunov scheme. 



I AUSMDV : S-K* ] 

pig, 7(b) A shock wave diffracting around a comer by AUS- 
MDV without entropy-fix. 


DENSITY 
Jmax- 71 
Kmax- 71 
CFL-1.00 
N - 124 
Tim#- 0.1307 



AUSMDV S-Ftx 


Fig. 7(c) A shock wave diffracting around a comer by AUS- 
MDV with entropy-fix. 
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CFL-1.00 N- 200 


4, NUMERICAL RESULTS AND DISCUS- 
SIONS 

Various numerical experiments, including a thermo- 
chemical nonequilibrium flow problem, are performed in 
order to verify the soundness of the AUSMDV with the 
entropy /shock-fix. 

4.1 One- Dimensional Shock- Tube Problems 

Probably, a theoretical drawback of the AUSMDV 
lies in its shock-capturing capability(SCC). The AUS- 
MDV is equipped with the properties: high-resolution 
for contact discontinuities; conservation of enthalpy for 
steady flows. They are obtained at the sacrifice of the 
theoretical basis of the SCC of Van Leer’s FVS, where 
the flux splittings are wisely designed to have degener- 
acy in subsonic region so as to build a stationary shock 
structure with no more than two interior zones. Hence 
the SCC of the AUSMDV should be assured by numer- 
ical experiments. For this purpose, several shock-tube 
problems were solved. 

Figure 8(a) is a solution by the AUSMDV scheme for 
a steady shock at Mach number of 25, where a primitive- 
variable-averaged intermediate state is inserted at the 
x = 0 point as an initial condition. This figure indicates 
the AUSMD holds almost the same SCC as that of Van 
Leer’s FVS, producing a monotonic and crisp shock pro- 
file. Figure 8(b) is the solution of the AUSM based on 
Mach number splitting[26], showing an overshoot. It is 
interesting, however, to see that the temperature has al- 
most no overshoot, suggesting that the overshoot of the 
AUSM is mainly due to density field(mass flux). Also it 
is noted that this overshoot depends on the initial con- 
dition. When the intermediate point was not inserted, 
the AUSM gave a smoother solution. 



Fig. 8(a) Shock-tube problem — stationary shock by AUS- 
MDV. 


AUSM[M-Sp«tting) 



Fig. 8(b) Shock-tube problem — stationary shock by 
AUSM. 

Another concern about SCC is a noise made by a 
slowly-moving shock wave. Roberts[18] proposed a test 
problem, in which a shock traverses at Mach number of 
3 in a computational domain of 100 grid points, passing 
a single cell in approximately 50 time steps at a CFL 
number of 0.95. Figures 9(a-f) show the strength of lin- 
ear and nonlinear waves after 2000 time steps by the 
AUSMDV, the AUSM, Van Leer’s FVS, Roe’s FDS, the 
HLLE and Osher’s scheme, respectively. Here, the Roe 
linearization was used to compute the wave strength. 
The Roe scheme and the HLLE, which belong to the 
Godunov-type scheme, produce a significant error be- 
hind the shock wave. Other schemes give much better 
solutions, although the AUSM has a small overshoot on 
ti — c wave. It may be noticed that the AUSMDV , AUSM 
and Van Leer’s FVS have undershoot on u wave, which 
was also a common problem for the Osher scheme dur- 
ing the computation. But this undershoot is confined 
near the shock wave, and is not considered serious. The 
Godunov- type schemes, on the other hand, leave strong 
perturbation with long wave-length behind the shock 
wave. 

The next problem is a moving contact discontinuity, 
of which condition is the same as in Section 3.3. Fig- 
ures 10(a-f) is the results by the AUSMDV, the AUSM, 
the Godunov scheme, the HLLE, Van Leer’s FVS, Van 
Leer’s FV/DS, respectively. The solution of the AUS- 
MDV is identical to that of the exact Riemann solver. 
Also Roe’s scheme and Osher’s scheme give the same re- 
sult. The solution of the AUSM is fairly good although 
slight oscillations are found in the pressure and veloc- 
ity fields. The other schemes are unacceptable. The 
HLLE is too dissipative, and the other FV/DS schemes 
diverged before 100 time steps, so Figures 10(e) and 
10(f) show the solutions at only 15-th time step. 
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X 


Fig. 9(a) Shock-tube problem — slowly moving shock wave 
by AUSMDV. 



X 


Fig. 9(b) Shock-tube problem — slowly moving shock wave 
by AUSM. 



x 


Fig. 9(c) Shock-tube problem — slowly moving shock wave 
by Van Leer’s FVS. 



x 


Fig. 9(d) Shock-tube problem — slowly moving shock wave 
by Roe scheme. 



x 


Fig. 9(e) Shock-tube problem — slowly moving shock wave 
by HLLE scheme. 



x 


Fig. 9(f) Shock-tube problem — slowly moving shock wave 
by Osher scheme. 
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Fig. 10(a) Shock-tube problem — moving contact disconti- Fig. 10(d) Shock-tube problem — moving contact disconti- 
nuity by AUSMDV. nuity by HLLE scheme. 




Fig. 10(b) Shock-tube problem — moving contact disconti- Fig. 10(e) Shock-tube problem — moving contact disconti- 
nuity by AUSM. nuity by Van Leer’s FVS. 




Fig. 10(c) Shock-tube problem — moving contact disconti- Fig. 10(f) Shock-tube problem moving contact disconti- 
nuity by Godunov scheme. nuity by Van Leer’s FV/DS. 
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Figures ll(a^c) show the solution of Sod’s problem: 
{p,P,v)l = (1> 1,0); (p,P,u)r = (l/ 8 » V 10 >°) by *ke 
AUSMDV, the AUSM and the Roe scheme, respectively. 
The AUSMDV gives almost the same result as the Roe 
scheme, while the AUSM shows slight overshoot and 
undershoot at the shock and the tail of expansion fan 
although it gives more accurate resolution of the fan. 

Finally, the robustness of the AUSMDV is veri- 
fied by the use of a shock-tube problem with strong 
expansion fans. The initial condition is (p,p, u)l = 
(1,2, — 2.5c&); ( p,p,u) R = (1, 0.5, 2.5c L ). Figures 12(a- 
c) are the solutions by the AUSMDV, the Godunov 
scheme and the HLLE, respectively. The AUSMDV 
gives a nice solution with no glitch. The Godunov 
scheme gives also a satisfactory solution, which was very 
similar to that by the AUSM, Van Leer’s FVS and Van 
Leer’s FV/DS, but there is a glitch at the sonic point. 
The HLLE scheme has no glitch, but its excessive dis- 
sipation gives a bump profile in the temperature. It is 
noted that Roe’s scheme and Osher’s scheme failed in 
this problem. 



Fig. 11(a) Sod’s shock- tube problem by AUSMDV. 



Fig. 11(b) Sod’s shock-tube problem by AUSM. 



Fig. 11(c) Sod’s shock- tube problem by Roe scheme. 



Fig. 12(a) Shock-tube problem — strong expansion by AUS- 
MDV. 



Fig. 12(b) Shock-tube problem — strong expansion by Go- 
dunov scheme. 
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Fig. 12(c) Shock-tube problem — strong expansion by HLLE 
scheme. 

4.2 Two-Dimensional Shear Flow 

The second problem is a two-dimensional shear flow 
through a constant-area duct using the Euler equa- 
tions. Two supersonic inlet flows divided at the cen- 
tral grid are assumed with the left boundary condition: 
(p,p, Af)upper = (1 j1i 2); (p»P> Slower = (10,1,1.1). 
Here a very coarse grid system of 10 x 10 was used. Fig- 
ures 13(a) shows the density contours and other quanti- 
ties solved by the AUSMDV, which maintains a strong 
shear layer. The Godunov scheme, Roe scheme, Os- 
her scheme and the AUSM also gave the same solu- 
tion. Figures 13(b-d) show the results by Van Leer’s 
FVS, Hand’s FV/DS, Van Leer’s FV/DS, respectivdy. 
These schemes can not keep the shear layer although 
Van Leer’s FV/DS shows a noticeable improvement over 
the Hanel’s FV/DS. It is noted that the HLLE scheme 
was as dissipative as the Van Leer FVS. 



p,u,p p,*,p P>">P 

AUSMDV : ApLf*d=PuK=JO CFL - 1 .00 N - 1000 


Fig. 13(a) Two-dimensional shear flow problem by AUS- 
MDV. 




p,*,p Pt**P 

Van Leer F.V.S. CFL - 1.00 N * 1000 


Fig. 13(b) Two-dimensional shear flow problem by Van 
Leer’s FVS. 
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p,*,p P**>P P>*>P 

Hanel89 F.V/D.S. CFL -1.00N - 1000 


Fig. 13(c) Two-dimensional shear flow problem by Hanel’s 
FV/DS. 




p,*,p p»*>p P’^P 

Van Laer/HaneOI CFL - 1 .00 N » 1 000 


Fig. 13(d) Two-dimensional shear flow problem by Van 
Leer’s FV/DS. 
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4.3 One-Dimensional Self- Similar Conical Flow 
Problem 

The third problem is also a standard problem: 1- 
D self- similar conical flow over a 10-degree half cone at 
hypersonic speed. This problem is used to check the ac- 
curacy of numerical schemes on the shock and boundary 
layer [36]. The flow conditions are Moo = 7.95; Re«> = 
4.2x10 s . Figures 14(a) and 14(b) show the profiles of 
pressure, temperature and tangential velocity solved by 
the AUSMDV and the Roe scheme, respectively. Since 
Pr = 1.0, the exact solution of adiabatic wall temper- 
ature is available and indicated by an arrow. The Roe 
scheme and the AUSMDV give similar results although 
the latter has a slightly steep curvature in the pressure 
field at the edge of the boundary layer. For this prob- 
lem, Van Leer’s FVS is known to thicken the boundary 
layer, as is shown in Fig. 14(c). Also Van Leer’s FV/DS 
gives a glitch[25] in Fig. 14(d). 


AUSMDV :ApuP*)=PiJC=10 &a»-°2 50 N * 5000 



Fig. 14(a) Conical flow problem by AUSMDV. 


Roe CFL-0.50 N- 5000 

E-Fix iwall-2 



Fig. 14(b) Conical flow problem by Roe scheme. 


Van Leer F.V.S. CFL-0.50 N- 5000 

IwalU 2 



Fig. 14(c) Conical flow problem by Van Leer’s FVS. 


Van Leer/Hanel91 CFL-0.50 N*= 5000 

iwalu 2 



Fig. 14(d) Conical flow problem by Van Leer’s FV/DS. 


AUSMDV :fj>^p*)=Pi*K=10 5TJrr°6 30 N - 10000 

E-Fix IwalU 2 



Fig. 14(e) Conical flow problem by second-order AUSMDV. 
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Roe CFL-0.30 N- 10000 

E-Rx »wa«- 2 



Van Leer/Hanel9l P F n*°o 30 10000 

fwaJI- 2 



Fig. 14 (g) Conical flow problem by second-order Van Leer 
FV/DS. 


The second-order results by the use of the MUSCL- 
type TVD scheme with minmod limiter, are shown in 
Figs. 14(e-g) for the AUSMDV, the Roe scheme and Van 
Leer’s FV/DS, respectively. The solutions of the AUS- 
MDV and the Roe scheme are almost identical while 
the kink of Van Leer’s FV/DS still remains although it 
is improved. 

4.4 Hypersonic Flow around a Blunt Body 

Inviscid and viscous hypersonic flows around a cylin- 
der were solved for the ideal gas. First the inviscid flow 
is presented at the free-stream Mach number of 6 with 
a 45x41 grid system and mirror wall condition being 
used. This grid system includes control volumes along 
the stagnation streamline. Figures 15(a-d) show the re- 
sults by the AUSMDV, the Godunov scheme, the Roe 
scheme and the Osher scheme, respectively. The density 
contours as well as the profiles of pressure, temperature 


and density along the stagnation stream line are shown. 
The AUSMDV clearly captures a strong shock wave with 
no overshoot. The Roe scheme shows a serious carbun- 
cle phenomenon, and the Osher scheme has a singularity 
in the density /temperature profile even though the pres- 
sure is correctly solved. It is noted that when the con- 
trol volumes were off stagnation stream line, the Osher 
scheme did not suffer from this phenomenon. The Go- 
dunov scheme is acceptable, however the density, which 
should increase by adiabatic compression along the stag- 
nation stream line, does not have a monotonic profile. 
The proposed shock-fix was applied to the Roe scheme, 
and the result is shown in Figure 15(e). The carbun- 
cle phenomenon is now cured, and the solution is now 
similar to that of the Godunov scheme. But the den- 
sity does not take the maximum value at the stagna- 
tion point. This situation is improved when a higher- 
order scheme is used. Figure 15(f) shows the result, in 
which the non-MUSCL fully-upwind second-order TVD 
scheme[5] was used. The non-conservation of enthalpy 
of the Roe scheme is compensated for by higher-order 
accuracy. On the other hand, the AUSMDV conserves 
the total enthalpy for the steady state, so it makes the 
stagnation temperature exact, and gives a mono tonic 
density profile. 

The full Navier-Stokes equations were solved using 
the same grid at Re = 1000, which results in a high cell 
Reynolds number of 25 at the wall. The wall tempera- 
ture was fixed at the free-stream temperature, and the 
no-slip boundary condition was used. The density on 
the wall was calculated by assuming no mass flux from 
the wall. Figures 16(a)-(d) show the solution of the the 
AUSMDV, Roe’s FDS, Van Leer’s FVS and Van Leer’s 
FV/DS, respectively. These figures show the contours 
of pressure and temperature, the profile of pressure and 
temperature along the stagnation stream line, and Stan- 
ton number on the cylinder surface. The solution by the 
AUSMDV in Fig. 16(a) is similar to that of the Roe 
scheme with shock-fix in Fig. 16(b). It is noticeable 
that the AUSMDV predicts a slightly higher Stanton 
number(heat flux) than that of the Roe scheme, indi- 
cating less amount of numerical dissipation in the AUS- 
MDV result. Figure 16(c) is a solution by Van Leer’s 
FVS, which produces a thick boundary layer, predicting 
much reduced heat flux. Figure 16(d) is a solution by 
Van Leer’s FV/DS, and exhibits a kink in the pressure 
field near the wall although the resolution of the bound- 
ary layer is improved. The same problem was solved on 
a finer grid(149 x 62, cell Reynolds number of 4 at the 
wall) by the use of the AUSMDV and the Roe scheme. 
The results are shown in Fig. 16(e) and 16(f), where 
both schemes predict almost the same heat flux. This 
grid dependency study indicates that the AUSMDV is 
capable of evaluating heat flux a s accurately as Roe’s 
approximate Riemann solver. 
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4.5 Hypersonic Flow with the Two- Temperature 
Model 

The AUSMDV was applied to a thermal-chemical 
nonequilibrium flow: a hypersonic flow around a cylin- 
der with diameter of lm. The cylinder is traveling at 
7km/sec(Mach 25.3) at an altitude of 85km. Park’s two- 
temperature model and Blottner’s 7-species chemical re- 
action model were used with a grid system of 32 x 50. 
Steady state solutions were obtained with the residual 
being reduced about three-order of magnitude by the 
point implicit time integration method. The in viscid 
flow solution is shown in Fig. 17, which includes the 
contours of pressure, translational-rotational tempera- 
ture and vibrational-electronic temperature along with 
profiles of these quantities as well as molar concentra- 
tions on the stagnation stream line. The AUSMDV can 
capture a strong shock wave crisply with no oscillation in 
the field of the pressure, translational/ vibrational tem- 
perature and chemical species concentrations. 



model around a cylinder by AUSMDV. 

The viscous flow was solved using the same flow 
condition and grid, although this grid is very coarse 
for viscous computation with the cell Reynolds num- 
ber near the wall being about 50. Figures 18(a), 18(b) 
and 18(c) are results by the AUSMDV, the generalized 
Roe scheme[37] with shock-fix and the Hanel FVS, re- 
spectively. The AUSMDV and the Roe scheme give 
similar solutions, while the Hanel FVS shows a very 
thick boundary layer due to the numerical dissipa- 
tion. The heat flux at the stagnation was 78.9KW/m 2 , 
76.2KW/m 3 , 13.6KW/m 3 for the AUSMDV, the Roe 
FDS and the Hanel FVS, respectively. This again con- 
firms that the AUSMDV is as accurate as the Roe 
scheme in evaluating heat flux. It is noted that, es- 
pecially in the nonequilibrium flow computations, the 
AUSMDV is much simpler and more efficient than the 
Roe-type scheme, because the latter needs the informa- 
tion of various physical quantities in order to construct 
a Roe-averaged Jacobian matrix[37, 38]. 



- 0.4 


Fig. 18 (a) Viscous hypersonic flow with the two- temperature 
model around a cylinder by AUSMDV. 



- 0.2 


Fig. 18 (b) Viscous hypersonic flow with the two-temperature 
model around a cylinder by Roe scheme. 



Fig. 18 (c) Viscous hypersonic flow with the two- temperature 
model around a cylinder by Hanel FVS. 
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5. CONCLUSION 

In this paper we have proposed a flux splitting 
scheme which is equipped with favorable properties: 
high-resolution for contact discontinuities; conservation 
of enthalpy for steady flows; numerical efficiency; appli- 
cability to chemically reacting flows. The results of var- 
ious numerical experiments including that of a thermo- 
chemical nonequilibrium flow indicate the soundness of 
the proposed scheme. In addition, a shock-fix is pre- 
sented, which cures the carbuncle phenomenon and is 
applicable to the Roe scheme as well as to the new 
scheme presented here. 
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